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HIGH RESOLUTION WIND 
MEASUREMENTS FOR OFFSHORE WIND 
ENERGY DEVELOPMENT 

CROSS-REFERENCE TO RELATED 5 

APPLICATIONS 

This application claims the benefit under 35 U.S.C. Section 
1 1 9(e) of the following co-pending and conmionly -as signed 
U.S. provisional patent application(s), which is/are incorpo- 
rated by reference herein: 

Provisional Application Ser. No. 61/173,916, filed on Apr. 

29, 2009, by Son Van Nghiem, and Gregory Neumann, 
entitled “HIGH-RESOLUTION WIND MEASUREMENTS x 5 
FOR OFFSHORE WIND ENERGY DEVELOPMENT,”. 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH AND DEVELOPMENT 

20 

The invention described herein was made in the perfor- 
mance of work under NASA contract number NAS7-1407, 
and is subject to the provisions of Public Law 96-517 (35 
USC 202) in which the Contractor has elected to retain title. 

25 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates generally to offshore wind 
energy development, and in particular, to a method, appara- 30 
tus, and article of manufacture for determining a fundamental 
dataset of offshore wind climatology at 1 -km posting, based 
on the decadal satellite scatterometer data record from the 
entire lifetime of the QuikSCAT mission. 

2. Description of the Related Art 35 

(Note: This application references a number of different 

publications as indicated throughout the specification by ref- 
erence numbers enclosed in brackets, e.g., [x]. A list of these 
different publications ordered according to these reference 
numbers can be found below in the section entitled “Refer- 40 
ences.” Each of these publications is incorporated by refer- 
ence herein.) 

Energy is a critical issue facing all from local to regional, 
national, and international levels. With recent development of 
wind tower technologies, offshore wind energy emerges as an 45 
important component of the global total renewable energy 
portfolio. Recognizing both the value and the complexity of 
offshore wind energy development, energy communities, 
commissions, and consortiums have been formed in many 
regions, which identify long-term high-resolution wind mea- 50 
surements as critical data to evaluate and select locations 
appropriate for wind farming. Thus, high-resolution wind 
measurements have a crucial role in addressing renewable 
energy and climate change. In this regard, to determine where 
to best locate an offshore wind farm, an accurate determina- 55 
tion/measurement of the location of offshore winds is neces- 
sary. 

The National Renewable Energy Laboratory (NREL) clas- 
sifies wind power density into seven classes [1] in which 
Class 4 or above is appropriate for wind farm development. 60 
For the estimation of the NREL power density classes, accu- 
rate long-term historical wind data are required. While model 
analyses are available in a number of selected regions, there 
are discrepancies among models and significant land con- 
tamination effects near land-sea boundaries [2]. Moreover, it 65 
is recommended that model estimates at any location be con- 
firmed by measurements. 


2 

One method used to confirm such measurements utilizes 
buoys at designated locations to obtained point measure- 
ments that may not represent wind conditions in the surround- 
ing area. In this regard, buoy measurements are obtained at a 
very limited number of locations and provide point data that 
do not represent the wind distribution over the region. For 
example, the buoy size itself is likely only 1 0-20 feet and only 
measures wind directly in the vicinity of the buoy. Accord- 
ingly, buoy calculations are not accurate or truly representa- 
tive of wind as the distance from the buoy increases. While 
additional buoys may be utilized, they cannot be placed every 
kilometer over the entire ocean surface and still only provide 
limited point data. If a buoy is deployed in every 1 km 2 across 
the world ocean, over 360 million buoys are needed. Assum- 
ing an impossibly low cost of $1 ,000 for each buoy per year, 
the total cost would amount to $3.6 trillion for a decade of 
data. Therefore, the extreme costs of manufacturing, deploy- 
ing, maintaining, and replenishing such a global massive 
buoy system over a decade are prohibitive let alone any other 
associated issues. 

Wind speed can be measured by satellite synthetic aperture 
radars; however, data are limited in accuracy, in time, and in 
space. Routine measurements of surface wind fields over 
global waters have been demonstrated by satellite scatterom- 
eter missions, where stable and accurate radars provide 
coarse resolutions [3-6] . In such missions, a wind scatterom- 
eter determines the normalized radar cross section (sigma) of 
the surface by transmitting a pulse of microwave energy 
towards the Earth’s surface and measuring the energy scat- 
tered backward (backscatter) to the scatterometer. By com- 
bining the sigma measurements from different azimuth 
angles, the near-surface wind vector over the ocean’s surface 
can be determined using a geophysical model function that 
relates wind and backscatter. However, due to the original 
coarse resolution (e.g., 25 km) and land contamination 
effects, standard wind products from the scatterometer mis- 
sions are not valid by 2 pixels or 50 km away from shore. 

While it is desirable for a wind farm to locate at some 
distance (-10 km) away from shore for higher and steadier 
wind speeds and because of issues related to ocean view, 
noise, navigation, security, and environmental impacts, wind 
tower technologies and bathymetry may limit the distance to 
within a few tens of km from shore (-50-100 km). Thus, the 
standard coarse low-resolution satellite wind products from 
satellite scatterometer data are neither useful nor applicable. 

To utilize satellite scatterometer data for offshore wind 
energy applications, the spatial resolution must be increased 
significantly. In this regard, attempts have been made to 
reconstruct high-resolution data based on low-resolution 
data. Algorithms for high-resolution data reconstruction from 
low-resolution data are traditionally based on the deconvolu- 
tion method, such as the Computer Tomography scan (CT 
scan) used in medical imaging, for which Hounsfield and 
Cormack earned the 1979 Nobel Prize in medicine. In this 
regard, the traditional deconvolution approach has been used 
to enhance resolution of radar data [7-9]. However, such 
approaches require that radar backscatter remains unchanged 
in each high-resolution pixel during the period of data acqui- 
sition used in the resolution enhancement process (similar to 
requiring a patient to stay still while the CT scan is made), and 
are not applicable to backscatter data acquired over ocean 
surface that are rapidly variable in time and in different azi- 
muth angles. In this regard, characteristics of water surface 
with wind-generated waves and roughness can change in a 
short time depending on wind variability. Moreover, the 
directional feature of waves causes differences in backscatter 
at different azimuth angles in data collected even within the 
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same orbit. These backscatter variations invalidate the funda- 
mental requirement of backscatter invariant for the deconvo- 
lution approach to be applicable. Therefore, a new method is 
necessary to account for backscatter change in the process to 
construct high-resolution data. 

In view of the above, what is needed is the capability to 
obtain high-resolution wind climatology for offshore wind 
energy applications (that account for backscatter change) and 
for use in other scientific research. 

SUMMARY OF THE INVENTION 

A new method to measure offshore/near-shore wind uti- 
lizes a unique dataset of offshore wind climatology at 1 -km 
posting from a decade of National Aeronautics and Space 
Administration (NASA) QuikSCAT satellite scatterometer 
data. High-resolution scatterometer results, as close as 10-15 
km to shore, are compared with model patterns of the 
National Renewable Laboratory (NREL) wind density 
classes to show potential areas for wind energy development 
in near-shore waters. High-resolution wind atlases can be 
made for various science studies and practical applications. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Referring now to the drawings in which like reference 
numbers represent corresponding parts throughout: 

FIG. 1 is an exemplary hardware and software environment 
used to implement one or more embodiments of the inven- 
tion; 

FIG. 2 illustrates a conical scanning configuration of the 
SeaWinds scatterometer aboard the QuikSCAT satellite; 

FIG. 3 presents the original Dense Sampling Method 
(DSM) backscatter image (inset) and a semi-transparent over- 
lay of the DSM image on a GOOGLE EARTH™ map of the 
island of Principe where clouds appeared as white areas in 
accordance with one or more embodiments of the invention; 

FIG. 4(a) illustrates a QuikSCAT backscatter signature 
pattern at 1-km posting in accordance with embodiments of 
the invention; 

FIG. 4(b) illustrates NREL power density estimated by the 
Hawaii Wind Working Group in accordance with one or more 
embodiments of the invention; 

FIG. 5 illustrates an example for the case of Hawaii to 
compare the high-resolution backscatter pattern with wind 
power density classes from model calculations in accordance 
with one or more embodiments of the invention; and 

FIG. 6 is a flow chart illustrating the logical flow for mea- 
suring wind in accordance with one or more embodiments of 
the invention. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

In the following description, reference is made to the 
accompanying drawings which form a part hereof, and which 
is shown, by way of illustration, several embodiments of the 
present invention. It is understood that other embodiments 
may be utilized and structural changes may be made without 
departing from the scope of the present invention. 

Overview 

One or more embodiments of the invention overcome the 
problems of the prior art by utilizes a transform operation that 
is applied on the mean and fluctuating parts of targeted pixels 
to construct data in a high-resolution grid at 1 -km postings for 
wind measurements over water surfaces such as oceans or 
lakes. 
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Hardware and Software Environment 

FIG. 1 is an exemplary hardware and software environment 
100 used to implement one or more embodiments of the 
invention. The hardware and software environment includes a 
5 computer 102 and may include peripherals. Computer 102 
may be a user/client computer, server computer, or may be a 
database computer. The computer 102 comprises a general 
purpose hardware processor 104 A and/or a special purpose 
hardware processor 104B (hereinafter alternatively collec- 
io tively referred to as processor 104) and a memory 106, suchas 
random access memory (RAM). The computer 102 may be 
coupled to other devices, including input/output (I/O) devices 
such as a keyboard 114, a cursor control device 116 (e.g., a 
mouse, a pointing device, pen and tablet, etc.) and a printer 
15 128. Additionally, satellite transmission data may be input 
into computer 102 directly or indirectly. In other words, radar 
based data (see detailed description below) may be received 
into computer 102 or received via another device and used by 
computer 102. 

20 In one embodiment, the computer 102 operates by the 
general purpose processor 104 A performing instructions 
defined by the computer program 110 under control of an 
operating system 108. The computer program 110 and/or the 
operating system 108 may be stored in the memory 106 and 
25 may interface with the user and/or other devices to accept 
input and commands and, based on such input and commands 
and the instructions defined by the computer program 110 and 
operating system 108 to provide output and results. 

Output/results may be presented on the display 122 or 
30 provided to another device for presentation or further pro- 
cessing or action. In one embodiment, the display 122 com- 
prises a liquid crystal display (LCD) having a plurality of 
separately addressable liquid crystals. Each liquid crystal of 
the display 122 changes to an opaque or translucent state to 
35 form a part of the image on the di splay in response to the data 
or information generated by the processor 1 04 from the appli- 
cation of the instructions of the computer program 110 and/or 
operating system 108 to the input and commands. The image 
may be provided through a graphical user interface (GUI) 
40 module 1 1 8 A. Although the GUI module 1 18 A is depicted as 
a separate module, the instructions performing the GUI func- 
tions can be resident or distributed in the operating system 
108, the computer program 110, or implemented with special 
purpose memory and processors. 

45 Some or all of the operations performed by the computer 
102 according to the computer program 110 instructions may 
be implemented in a special purpose processor 104B. In this 
embodiment, the some or all of the computer program 110 
instructions may be implemented via firmware instructions 
50 stored in a read only memory (ROM), a programmable read 
only memory (PROM) or flash memory within the special 
purpose processor 104B or in memory 106. The special pur- 
pose processor 104B may also be hardwired through circuit 
design to perform some or all of the operations to implement 
55 the present invention. Further, the special purpose processor 
104B may be a hybrid processor, which includes dedicated 
circuitry for performing a subset of functions, and other cir- 
cuits for performing more general functions such as respond- 
ing to computer program instructions. In one embodiment, 
60 the special purpose processor is an application specific inte- 
grated circuit (ASIC). 

The computer 102 may also implement a compiler 112 
which allows an application program 110 written in a pro- 
gramming language such as COBOL, Pascal, C++, FOR- 
65 TRAN, or other language to be translated into processor 104 
readable code. After completion, the application or computer 
program 110 accesses and manipulates data accepted from 
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I/O devices and stored in the memory 106 of the computer 
102 using the relationship s and logic that was generated using 
the compiler 112. 

The computer 102 also optionally comprises an external 
communication device such as a modem, satellite link, Eth- 5 
emet card, or other device for accepting input from and pro- 
viding output to other computers. 

In one embodiment, instructions implementing the operat- 
ing system 108, the computer program 110, and the compiler 
112 are tangibly embodied in a computer-readable medium, 10 
e.g., data storage device 120, which could include one or 
more fixed or removable data storage devices, such as a zip 
drive, floppy disc drive 124, hard drive, CD-ROM drive, tape 
drive, etc. Further, the operating system 108 and the computer 
program 110 are compri sed of computer program instructions 1 5 
which, when accessed, read and executed by the computer 
102, causes the computer 102 to perform the steps necessary 
to implement and/or use the present invention or to load the 
program of instructions into a memory, thus creating a special 
purpose data structure causing the computer to operate as a 20 
specially programmed computer executing the method steps 
described herein. Computer program 110 and/or operating 
instructions may also be tangibly embodied in memory 106 
and/or data communications devices 130, thereby making a 
computer program product or article of manufacture accord- 25 
ing to the invention. As such, the terms “article of manufac- 
ture,” “program storage device” and “computer program 
product” as used herein are intended to encompass a com- 
puter program accessible from any computer readable device 
or media. 30 

Of course, those skilled in the art will recognize that any 
combination of the above components, or any number of 
different components, peripherals, and other devices, may be 
used with the computer 102. 

Although the term “user computer” or “client computer” is 35 
referred to herein, it is understood that a user computer 102 
may include portable devices such as cell phones, notebook 
computers, pocket computers, or any other device with suit- 
able processing, communication, and input/output capability. 
Software Embodiments 40 

Embodiments of the invention may be performed using the 
computer hardware and software environment described 
above. More specifically, the computer program 110 may 
implement and provide the ability to construct data in a high- 
resolution grid based on a transform applied to various data 45 
while accounting for backscatter change. Initially, one begins 
with low resolution data available from a satellite scatterom- 
eter. Thereafter, a dense sampling method (DSM) is used to 
construct the data in a high-resolution grid. 

Satellite Scatterometer Data 50 

Among several global scatterometer datasets, the SeaW- 
inds scatterometer mission with the QuikSCAT (QS) satellite 
has been providing the most stable and accurate backscatter 
data across the world on a daily basis for more than a decade. 

QS was launched from Vandenberg Air Force Base in Cali- 55 
fornia in June 1999, and continued to operate in a low-earth 
orbit along sun- synchronous passes until November 2009. 

QS acquired data with a scanning pencil -beam antenna across 
large swaths over the Earth’s surface [10]. QS collected glo- 
bal data covering 90% of the Earth’s surface each day [11]. 60 
The QS footprint (two-way, half-power, full beam width) has 
an elongated or ‘egg’ shape (FIG. 2) of about 25 km in 
azimuth by 37 km in range. In this regard, FIG. 2 illustrates a 
conical scanning configuration of the Sea Winds scatterom- 
eter. The outer scan has a swath width of 1 800 km for the 65 
vertical polarization and the inner swath is 1400 km for the 
horizontal polarization, which are centered along the nadir 


track. Moreover, the transmitted radar pulse is modulated 
(chirped), and the received pulse is Doppler compensated 
before it is processed through a fast-Fourier-transform stage 
to achieve a sub-footprint resolution. This linear frequency 
modulation chirp function is capable of generating sub -foot- 
prints or thin “slices” of 25 km in azimuth with a commend- 
able range resolution of 2- 1 0 km with a nominal value at 6 km 
[11]. This QS dataset is used as the low -resolution data in 
embodiments of the invention and is input into computer 
system 102. 

Dense Sampling Method 

The new method, called the Dense Sampling Method 
(DSM), is based on the Rosette Transform applied on the 
mean and fluctuating parts of each targeted pixel to construct 
the data in a high-resolution grid at 1-km posting. DSM has 
been successfully applied to observed urban and suburban 
environments [12]. Unlike the deconvolution method, DSM 
allows the use of the low-resolution data collected by the 
satellite scatterometer from space to construct high resolution 
data over targets with varying signatures such as on water 
surfaces roughened by winds. 

Consider an x-y coordinate surface where the center of a 
radar footprint is located at (x^,yj. In the formulation below, 
( x £,>T) i s se t a * (0,0) for simplicity without loss of generality, 
and all backscatter quantities are in the linear domain (as 
opposed to the decibel or dB domain that is logarithmic and 
non linear). The radar gain pattern is approximately modeled 
with an elongated Gaussian beam tilted at an azimuth angle of 
(|) z as defined by the following expression 


G(<pi, x, y) = exp - 


x cos (pi + y sin0,- 


M 


-x sin (pi + y cos(p l 


41 


(1) 


where parameters a and b are determined by the full beam- 
width W x in azimuth and W v in range, respectively with the 
equations below 


w x w y (2) 

2V^2 2sfhi2 


At the azimuth angle of (|) z =0, the Gaussian beam has the 
simple form of 

G(<^ i ,x,y)=exp(-x 2 /a 2 -y 2 /b 2 ) (3) 

Backscatter cr 0 ((|) z , t z , x, y) at location (x,y), measured at 
time t z in an azimuth direction of <|> z ., is derived from the 
antenna power beam G((|) z , x, y), transmitted power P rz , 
received power P^ z , wavelength X, and distance R between the 
satellite and the targeted area A far away from the scatterom- 
eter by the radar equation 


r r 

PRi = (4tt) 3 p 4 J J A dxd y G ^i' x ' ywti* y) 


(4) 


Because the footprint covers many locations with different 
backscatter values in a targeted area, a single radar measure- 
ment at azimuth angle ^ can be characterized with an average 
backscatter value a 0 (( |) z , t z ) to allow it to be factored out of the 
integration so that 
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P Ti A 2 
(4 n?& 


^o(0i 


-«/X 


dxdyG{(p jv, y) 


( 5 ) 


8 

ntij 

ij/x 


-continued 

N 

dxdyG((p x, ;y)£(0 r,-, a:, y) 


In this formulation, ct 0 (())„ t f ) represents the measured back- Each of the two terms in (1 0) can be examined to describe 

scatter and other backscatter values are calculated. Substitut- i ts me aning For the first term, denoted as a 0M , the order of the 

ing (5) into (4), the averaged backscatter can be expressed as: summation and the double integration are interchanged, and 

the mean part cr 0 (x, y) is factored out of the summation since 
it is independent of angle <|> z . at time t z so that becomes: 


o-o((pi, ti ) = f dxdyG(<p;, x, y)cr 0 (<p t h x, y) 

where the integrated gain-area product has been defined as 15 

r A = JJ dxdyG(<p;,x, y) 


r.IL‘ 


(ID 


\<r o(x, y) 


dxdyG(x, y)cr 0 (x, y) 


( 7 ) 


20 In ( 1 1 ), the equivalent composite power gain is defined as : 


The formulation given by (6) is similar to the two-dimen- 
sional Radon transform used in tomography imaging when 
G((|)., x, y) is replaced with the Dirac delta function [13]. 
However, as described above, such a radon transform is used 25 
when the backscatter remains unchanged and cannot be used 
for targets with varying signatures such as on water surfaces 
roughened by winds. 

An ensemble average of N samples of backscatter data may 30 
be measured at azimuth angle (|) z at time t z densely populated 
within area A. From many measurements with their footprint 
centers collocated at (x^,y t )=(0,0), the backscatter ensemble 
a 0 is represented by the following equation 


k§ji 


dxdyG{<p x, y)cr 0 ((p r ; , x, y) 


( 8 ) 


Similar backscatter ensembles are carried out at all other 
locations (x^y*), densely posted in each 30 arcsec in latitude 
and in longitude (about 1 km near the equator). As such, DSM 
includes data densely sampled in terms of both azimuth diver- 
sity and posting locations. Note that data with their centers 45 
outside of each posting are not included at all in the ensemble 
at that posting. This is equivalent to imposing a binary 
weighting function W(x^,y^) with a value of 1 inside and 0 
outside the 30-arcsec posting area centered at (x^y^). 

The target may be decomposed into two different parts: (1 ) 50 
the average represented as a 0 (x, y); and (2) the changing or 
fluctuation represented by €((|) z , t f , x, y). In this regard, at each 
coordinate (x,y), backscatter ct 0 ((|) z , t., x, y) at the posting 
(x*,y*M0,0) is composed of a mean part of cr 0 (x,y), averaged 
over all (|) z and all t z , and a zero-mean fluctuation part of €((|) z , 
t z ., x, y) variable in (|). at different time t z ., so that 

tf o(<h> tp x,y ) = d 0 (x,y )+£ (<J)j, tpx,y) (9) 

Using (9) in (8) and separating the mean and fluctuation parts, 60 
one can obtain: 


N 

G(x, y) = £ 
1 = 1 


G(<p h x, y) 
N 


(12) 


As expressed in (12), the composite gain is a linear com- 
bination of the elongated patterns of the thin slices at different 
azimuth angles. Accordingly, from the decomposed target, 
the measurements can be reconstructed at a higher resolution. 
For a figurative description of (12), consider a simple case 
where G(x,y) consists of two perpendicular slices with their 
center collocated at (x Ar ,y Ar )=(0,0). In this case, the elongated 
patterns of the two thin slices form a composite pattern simi- 
lar to a four-petal rosette. The four-petal pattern has a high 
gain around the center and the gain gradually tapers into the 
four-petal sides. When G(x,y) includes multiple thin slices 
over different azimuth angles, the composite pattern becomes 
similar to a multi-petal rosette. As such, the transformation 
defined by (1 1) and (12) is called the Rosette Transform. For 
a large number of collocated thin slices collected in the full 
range of azimuth angles, the composite gain G(x,y) in the 
Rosette Transform approaches an isotropic pattern with the 
gain uniformly distributed around the peak at the center. 

Since all backscatter samples are overlapped near the cen- 
ter of each posting and are separated at the far ends of the 
individual elongated power pattern, the composite gain is 
preserved at the center while it is more weighted down farther 
away from the center. In this method, the multi-azimuth com- 
position of the thin slices forms the equivalent gain G(x,y) 
that achieves a higher resolution approximately at the scale of 
the smallest dimension of the thin slice footprint for the mean 
backscatter at each posting (every 30 arcsec in latitude and in 
longitude). 

The second term in (12), represented by , can be 

rearranged in the form of 



fit’ll 


G((pj, x, y) 
N 


x, y) 


(13) 


o-o = j f dxd yG((Pi , y)o-o(x, y) H 


( 10 ) 

65 where the fluctuation term €((|) z , t z , x, y) cannot be factored out 
of the summation because both the fluctuation term and the 
individual gain are azimuthally dependent. 
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It is noted that, while e((|) z , t f , x, y) has a zero mean as 
defined by (9), the ensemble average in (13) of g(c|> z . ? t t , x, y), 
which is individually weighted by a different gain factor G((|) z , 
x, y) at each azimuth angle (|);, is not necessarily zero in 
general. 

For a small number of samples, the mean part is not repre- 
sentative and the residual makes the composite result of back- 
scatter noisy. By a sacrifice of the temporal resolution of 
scatterometer data over a long period of time, such as over a 
year or over multiple years or over the QS lifetime, one can 
obtain a stable mean of backscatter at each posting, and the 
residual becomes small as it can be positive at one azimuth 
angle and negative at another azimuth angle. Moreover, (13) 
indicates that a fluctuation outlier is suppressed directly by 
the sample number N. This is because N appears in the 
denominator inside the summation in (13) for each fluctua- 
tion term €(c|) z , t z , x, y) corresponding to index i. Thereby, an 
outlier fluctuation term due to an infrequent abnormal case is 
divided by N when it is averaged together with many other 
samples at different azimuth angles. Since a long time period 
is needed to obtain sufficient data for a stable mean with a 
small residual, a necessary condition of DSM is that the 
targeted area maintains persistent characteristics or climato- 
logical patterns during each long period of the ensemble 
averaging process. 

Note that (10) allows variability in azimuth and in time 
within each targeted area, such as the azimuth asymmetry of 
wind-driven waves on ocean surface. Natural occurrences, 
including wind, snow, rains, floods, droughts, or diurnal tem- 
perature changes, affect the geometric and electromagnetic 
properties of the environments causing changes in radar back- 
scatter. These changes are included in the residual term given 
by (13). The variability of each backscatter ensemble can be 
characterized by its variance or its standard deviation to rep- 
resent how much change occurred in each posting. The stan- 
dard deviation is related to variations (such as wind speed 
distribution) and can be used to characterize the variability. 
Furthermore, the general DSM formulation allows the use of 
data at both vertical (W) and horizontal (HH) polarizations, 
and their different signature is included in the standard devia- 
tion. 

The new concept of DSM presented above is different from 
the deconvolution approach traditionally used to enhance 
resolution of radar data ([7], [8], [18], and [19]). DSM is 
simply based on a linear composition of a set of thin-slice 
beam measurements as indicated by [12] to obtain a sharper 
and higher-resolution result by the multi-azimuth composi- 
tion instead of deconvolving overlapping beams. DSM sup- 
presses fluctuations in backscatter of a stationary target due to 
azimuth asymmetry and environmental changes, while such 
changes together with noise cause spurious artifacts ampli- 
fied by the non-linear deconvolution process. 

In the deconvolution approach, [8] assume azimuthally 
independent radar echo, which is applicable to an isotropic 
area such as a forested area with vegetation randomly ori- 
ented in azimuth, and such assumption is not valid for an 
azimuthally asymmetric area (such as ocean surface). Espe- 
cially for wind over an ocean or a lake, backscatter can be 
significantly different even within one single orbit (fore and 
aft looks) dues to the azimuth asymmetry. The deconvolution 
approach further requires that backscatter of the targeted area 
remains unchanged during the time period of data acquisition 
for the resolution enhancement ([9]). Such requirements and 
the azimuth asymmetry invalidate the resolution enhance- 
ment result by the deconvolution method when a transient 
effect (e.g., wind, snowmelt, precipitation, flood, diurnal 
effects, etc.) causes a backscatter change in the targeted area. 


Results and Verification 

A prototype of the above method has been developed to 
obtain initial high-resolution backscatter to validate the appli- 
cability of the method to QS data and to verily the accuracy of 
5 the new results for high-resolution backscatter. There are 
three aspects to confirm: (1) can the geographic target be 
correctly located?; (2) can the shape of the target be recon- 
structed accurately?; and (3) can the surface area be calcu- 
lated accurately? 

to For the verification purpose, a target can be selected with 
backscatter higher than its surrounding environment, such as 
a small island isolated in an ocean where geophysical condi- 
tions are stable and no significant change of the coastline has 
occurred since the launch of QS in 1 999. A candidate for this 
1 5 test is the Principe Island, located around 1 . 62 ° N and 7 .40° E 
off the central west coast of Africa. This island is famous as 
the site where Eddington experimentally proved Einstein’s 
theory of relativity using a solar eclipse in 1919 [14]. 

FIG. 3 presents the original DSM backscatter image (inset) 
20 and a semi-transparent overlay of the DSM image on a 
GOOGLE EARTH™ map of the island of Principe where 
clouds appeared as white areas. The DSM backscatter thresh- 
old of -12.5 dB (black contour in the inset, and white contour 
in the map overlay) is consistent with the island extent. 

25 The island has backscatter higher than that of the surround- 

ing ocean as seen in the scale from -10 dB (red 302 ) on the 
island down to -15 dB at the dark edge 304 in the ocean 
around the island. Note that the coast line contour at - 1 2 .5 dB 
is simply the average between -10 dB of the island and -15 
30 dB of the ocean. The 30-arcsec posting gives a smooth tran- 
sition between two adjacent areas with low and high back- 
scatter values, enabling a consistent identification of the 
boundary between the two areas. The result in FIG. 3 indi- 
cates that the backscatter threshold of -12.5 dB corresponds 
35 with the extent of the island within a couple of km. Within this 
threshold, the area calculated from the DSM image is 130 
km 2 , which is very close to a true area of 136 km 2 of the 
Principe Island. The general shape of the DSM image is also 
consistent with the shape of the island, which is elongated in 
40 the north-south direction. Furthermore, the geolocation of the 
DSM image coincides well with the map of the island. This 
case of a real target verifies the capability of DSM in delin- 
eating the high-resolution extent of a target and the accuracy 
of the DSM geolocation derived from QS LIB data. In another 
45 case [15], it can be shown that DSM is able to identify and 
accurately locate the Nukunonu Atoll in the South Pacific 
Ocean, a circular strip of land with a width of about 200 m 
located approximately 1280 km to the east-northeast of Fiji. 
Procedure for Wind Energy Estimate 
50 The procedure for offshore wind energy applications 
includes: (a) high-resolution backscatter data processing, (b) 
calculations of wind fields from QS backscatter empirical 
formulation, (c) verification of wind results with available 
surface measurements from buoys and ships, and (d) wind 
55 power estimates using the satellite high-resolution wind cli- 
matology dataset. The implementation of the high-resolution 
algorithm needs to include data processing for long-term data 
(~1 0 years) and for data partition in time to determine diurnal, 
seasonal, and interannual variabilities. For item (b), empirical 
60 geophysical model functions [5, 11, and 16] can be used to 
calculate high-resolution wind speed from high-resolution 
backscatter obtained in step (a). 

FIG. 4(a) illustrates a QS backscatter signature pattern at 
1 -km posting in accordance with embodiments of the inven- 
65 tion, and FIG. 4(b) illustrates NREL power density estimated 
by the Hawaii Wind Working Group [17]. Similarly, FIG. 5 
illustrates an example for the case of Hawaii to compare the 
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high-resolution backscatter pattern with wind power density 
from model calculations [17] in accordance with one or more 
embodiments of the invention. FIG. 4 (a) and the upper panel 
of FIG. 5 are for a QS backscatter (left scale in dB) signature 
pattern at 1 -km posting derived from DSM. FIG. 4(b) and the 5 
lower panel of FIG. 5 are for power density estimated by the 
Hawaii Wind Working Group [17] with the NREL classes 
marked in FIG. 5 by the white numbers. Ocean bathymetry is 
represented by the shaded relief in FIG. 4(b) and the lower 
panel of FIG. 5. 10 

The patterns of DSM backscatter and of model-calculated 
NREL power density classes compare consistently overall. 
DSM backscatter has land effects within 10-15 km from 
shore, where DSM results should not be used and such water 15 
area adjacent to the coastline will not be useful for wind 
energy development anyway as pointed out earlier. The com- 
parison reveals discrepancies such as anomalous lobes of 
high wind energy in the model estimates that are not consis- 
tent with the DSM pattern from QS satellite data. The results 20 
illustrate the deficiency in model results, containing biases 
and artifacts so that the Hawaii Wind Working Group itself 
recommends that actual measurements should be used to 
obtain an accurate estimate of wind energy. 

Examining FIG. 5, one may note that in the upper panel, 25 
between Hawaii and Maui, the DSM pattern protrudes 
approximately toward the southwest direction without a back 
lobe of high wind. However, the corresponding area in the 
lower panel indicates a different wind energy density pattern 
with a back lobe of high energy (grey area of class 4 to the 30 
right of class 6 between Hawaii and Maui) that is inaccurate. 
Instead, the patterns determined using the methods of the 
invention illustrated in the upper panel more accurately 
reflect the true wind patterns around the Hawaiian Islands. 35 
Further, examining FIG. 4(b) from the model estimate, one 
may note the dark shaded lobe around the southern tip of 
Hawaii indicating high wind energy which is not depicted in 
FIG. 4(a) from the DSM satellite measurement. 

In addition to the accuracy described above, it may further 40 
be noted that if one were to utilize prior art techniques, Maui 
itself is merely a couple footprints and such result are there- 
fore neither useful nor applicable. 

Logical Flow 

FIG. 6 is a flow chart illustrating the logical flow for mea- 45 
suring wind in accordance with one or more embodiments of 
the invention. 

At step 600, data at a first resolution (i.e., low resolution 
data) that has been collected by a satellite scatterometer in 
space is obtained. In this regard, the data may be obtained 50 
from the QuikSCAT scatterometer mission over a ten-year 
period of time (or other active or passive radiation data). The 
low resolution data may be represented by a radar equation: 

55 

Pm = f dxdyGtyi, x, y)o- 0 ((f>i , r f , y) 

wherein, cr 0 ((|) z ., t z , x, y) comprises backscatter at location (x,y) 60 
measured at time t z in an azimuth direction of <|> z -, G((|)., x, y) 
comprises an antenna power beam, comprises transmitted 
power, P Ri comprises received power, X comprises wave- 
length, and R comprises a distance between the satellite and a 
targeted area A. 65 

At step 602, one or more thin slices of the data are deter- 
mined. 


At step 604, a collocation of the data slices at a grid cell 
center are determined. Step 604 results in one or more 
ensembles of collocated data slices. 

At step 606, the collocated data slices are decomposed into 
a mean/average part and a fluctuating part. The decomposing 
of the collocated data slices is represented by: 

c* 0 (<tU t i ,x,y)=o 0 {x,y)+^ i , t p x,y) , 

wherein at each coordinate (x,y), backscatter a 0 (c|) z , t z , x, y) at 
a posting (x jt ,y jt )=(0,0) is composed of the mean part a 0 (x, y), 
averaged over all (|) z and all t z , and a zero -mean fluctuation part 
of €((|) z , t z , x, y) variable in at different time t z . 

At step 608, the data is reconstructed into a second resolu- 
tion (i.e., high resolution data) including the mean part and a 
residue of the fluctuating part. The reconstructing of the high- 
resolution result of the mean part, using a rosette transform, is 
defined by: 


^0 M 


dxdyG(x, y)v 0 (x, y) 


F o(x, y) 


wherein a QM represents the mean value at the high resolution, 
T a defines an integrated gain-area product, G(<|) z ., x, y) defines 
an antenna beam at an azimuth angle ^ at point (x,y), and an 
equivalent composite power gain is defined as: 


N 

G(x, y) = J] 

i = 1 


G((Pi,x, y) 
N 


At step 610, a wind measurement (i.e., high resolution 
wind results) is determined from the high-resolution data 
construction using a wind model function. 

At step 610, a description of the high-resolution wind 
results are output (e.g., to a graphical user interface, file, or 
other mechanism configured to use/display the results). 
Conclusion 

This concludes the description of the preferred embodi- 
ment of the invention. The following describes some alterna- 
tive embodiments for accomplishing the present invention. 
For example, any type of computer, such as a mainframe, 
minicomputer, or personal computer, or computer configura- 
tion, such as a timesharing mainframe, local area network, or 
standalone personal computer, could be used with the present 
invention. 

The foregoing description of the preferred embodiment of 
the invention has been presented for the purposes of illustra- 
tion and description. It is not intended to be exhaustive or to 
limit the invention to the precise form disclosed. Many modi- 
fications and variations are possible in light of the above 
teaching. It is intended that the scope of the invention be 
limited not by this detailed description, but rather by the 
claims appended hereto. 
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What is claimed is: 

10 1. A computer implemented method for measuring high- 

resolution wind, comprising: 

obtaining, in a computer, data at a first resolution collected 
by a satellite scatterometer, wherein the data comprises 
spatially dependent backscatter that is variable in time 
15 and in different azimuth angles; 

determining, in the computer, one or more thin slices of the 
data that are variable in time and in different azimuth 
angles; 

determining, in the computer, a collocation of the one or 
20 more data slices , that are variab le in time and in different 

azimuth angles, at each of one or more grid cell centers 
to obtain one or more ensembles of collocated data 
slices; 

decomposing, in the computer, the collocated data slices 
25 into a mean part, and a fluctuating part that is variable in 

time and in different azimuth angles; 
reconstructing, in the computer, the data at a second reso- 
lution from the mean part and a residue of the fluctuating 
part; 

30 determining, in the computer, a wind measurement from 
the data at the second resolution using a wind model 
function; and 

outputting, from the computer, a description of the wind 
measurement. 

35 2. The computer implemented method of claim 1, wherein 

the data at the first resolution is collected from a QuikSCAT 
scatterometer mission over a ten-year period of time. 

3. The computer implemented method of claim 1, wherein 
the data at the first resolution is represented by a radar equa- 
40 tion: 


p Ti x 2 r r 

PRi = (4 7T) 3 /? 4 J J A dxdyG ^ h X ’ >) 

45 

wherein, cr 0 ((|) z , t z , x, y) comprises spatially dependent 
backscatter, that is variable in time and in different azi- 
muth angles, to be decomposed into the mean part and 

50 the fluctuating part at location (x, y) measured at time t z 
in an azimuth direction of <(>,-, G(c() z , x, y) comprises an 
antenna power beam, P rz comprises transmitted power, 
P^. comprises received power, X comprises wavelength, 
and R comprises a distance between the satellite and a 

55 targeted area A. 

4. The computer implemented method of claim 1, wherein 
the decomposing of the collocated data slices is represented 

by: 

60 <*o(<L h x, y)=o 0 (x, D+e(<L h X, y), 

wherein at each coordinate (x, y), backscatter cr 0 ((|) z ., t., x, y) at 
a posting (x*, yjt) = (0,0) is composed of the mean part a 0 (x, y), 
averaged over all (|) z and all t ., and a zero-mean fluctuation part 
of g( 4> z , t z , x, y) variable in (|) z at different time t z . 

65 5 . The computer implemented method of claim 1, wherein 

the reconstructing the data at the second resolution from the 
mean part utilizes a rosette transform that is defined by: 
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dxdy\ 


g(0«, y) 

i N 


c 0 (x, y ) 


dxdyG(x , ^)ct 0 (j:, y) 


wherein a 0A ^ represents the mean part at the second resolution 
redefines an integrated gain-area product, G(<|), x, y) defines 
a Gaussian beam at an azimuth angle (|) z at point (x, y), and an 
equivalent composite power gain is defined as: 


N 

G(x,y) = Yj 

i= 1 


G{(pj,x, y) 
N 


6. An apparatus for measuring wind in a computer system 
comprising: 

(a) a computer having a memory; 

(b) an application executing on the computer, wherein the 

application is configured to: 

(i) obtain data at a first resolution collected by a satellite 
scatterometer wherein the data comprises spatially 
dependent backscatter that is variable in time and in 
different azimuth angles; 

(ii) determine one or more thin slices of the data that are 
variable in time and in different azimuth angles; 

determining a collocation of the one or more data slices, 
that are variable in time and in different azimuth 
angles, at each of one or more grid cell centers to 
obtain one or more ensembles of collocated data 
slices; 

(iii) decompose the collocated data slices into a mean 
part and a fluctuating part; 

(iv) reconstruct the data at a second resolution from the 
mean part and a residue of the fluctuating part; 

(v) determine a wind measurement from the data at the 
second resolution using a wind model function; and 

(vi) output a description of the wind measurement. 

7. The apparatus of claim 6, wherein the data at the first 
resolution is collected from a QuikSCAT scatterometer mis- 
sion over a ten-year period of time. 

8. The apparatus of claim 6, wherein the data at the first 
resolution is represented by a radar equation: 


Pri = f f dxdyGtyi, x, y)o- 0 (<t>i, t h x, y) 

wherein, cr 0 ((|) z , t z , x, y) comprises spatially dependent back- 
scatter that is variable in time and in different azimuth angles, 
to be decomposed into the mean part and the fluctuating part 
at location (x, y) measured at time t z in an azimuth direction of 
(|) z , G(j) z ., x, y) comprises an antenna power beam, P rz . comprises 
transmitted power, p» comprises received power, X com- 
prises wavelength, and R comprises a distance between the 
satellite and a targeted area A. 

9. The apparatus of claim 6, wherein the application is 
configured to decompose the collocated data slices based on 
an equation: 

a 0 (<|u t u x, y)=o 0 {x, y)+e(<l>* ti, x, y), 

wherein at each coordinate (x, y), backscatter cr 0 (<|) z -, t z -, x, y) at 
a posting (x*, y*)=(0,0) is composed of the mean part a 0 (x, y) 
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averaged over all <|) z and all t z , and a zero -mean fluctuation part 
of €(<|) z ., t z ., x, y) variable in <|) z . at different time t z .. 

10. The apparatus of claim 6, wherein the application is 
configured to reconstruct the data at the second resolution 
from the mean part utilizing a rosette transform that is defined 
by 


10 


°~0M '■ 


7/Hg 


G((pj, x , y) 
N 


ko(*. y) 


dxd yG(x, y)cr o(x, y) 
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wherein a 0M represents the mean part at the second resolution 
T a defines an integrated gain-area product, G((|), x, y) defines 
a Gaussian beam at an azimuth angle <^ i at point (x, y), and an 
equivalent composite power gain is defined as: 

20 


N 

G(x, y) = Y j 

i = 1 

25 


G(<Pi,x, y) 
N 


11. A non-transitory computer readable storage medium 
encoded with computer program instructions which when 
accessed by a computer cause the computer to load the pro- 

30 gram instructions to a memory therein creating a special 
purpose data structure causing the computer to operate as a 
specially programmed computer, executing a method of mea- 
suring wind, comprising: 

(a) obtaining, in the specially programmed computer, data 
at a first resolution collected by a satellite scatterometer 
wherein the data comprises backscatter that is variable in 
time and in different azimuth angles; 

(b) determining, in the specially programmed computer, 

40 one or more thin slices of the data that are variable in 

time and in different azimuth angles; 

(c) determining, in the specially programmed computer, a 
collocation of the one or more data slices, that are vari- 
able in time and in different azimuth angles, at each of 

45 one or more grid cell centers to obtain one or more 
ensembles of collocated data slices; 

(d) decomposing, in the specially programmed computer, 
the collocated data slices into a mean part and a fluctu- 
ating part that is variable in time and in different azimuth 
angles; 

(e) reconstructing, in the specially programmed computer, 
the data at a second resolution from the mean part and a 
residue of the fluctuating part; 

55 (f) determining, in the specially programmed computer, a 

wind measurement from the data at the second resolu- 
tion using a wind model function; and 

(g) outputting, in the specially programmed computer, a 
description of the wind measurement. 

12. The non-transitory computer readable storage medium 
of claim 11, wherein the data at the first resolution is collected 
from a QuikSCAT scatterometer mission over a ten-year 
period of time. 

65 13. The non-transitory computer readable storage medium 

of claim 11, wherein the data at the first resolution is repre- 
sented by a radar equation: 
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Pm = 


Pti£ 

(4tt )-’/<" 


If, 


dxdyG((pi, 


x, y)o- 0 (<f>i, t h x, y) 


15. The non-transitory computer readable storage medium 
of claim 11, wherein reconstructing the data at the second 
resolution from the mean part utilizes a rosette transform that 
is defined by: 


5 

wherein, cr 0 (4> z , t z , x, y) comprises spatially dependent back- 
scatter that is variable in time and in different azimuth angles, 
to be decomposed into the mean part and the fluctuating part 
at location (x, y) measured at time t z in an azimuth direction of 
(|) z , G((|) z , x, y) comprises an antenna power beam, com- 1C 
prises transmitted power, p* comprises received power, A, 
comprises wavelength, and R comprises a distance between 

the satellite and a targeted area A. where ‘ n ^presents the mean part at the second resolu- 

tion, 1 A defines an integrated gain-area product, G(q) ., x, y) 
14. The non-transitory computer readable storage medium 15 defines a Gaussian beam at an azimuth angle (|>, at point (x, y), 
of claim 11, wherein the decomposing ot the collocated data and an equivalent composite power gain is defined as: 
slices is represented by: 


°~0M '■ 


r.ff. 

rJL 


dxdyl 


i G(0/, y) 
i N 


k o(x, y) 


dxdyG(x, y)(r 0 (x , y) 


cr 0 (<tu t h x, y)=o 0 (x, v)+e(<k, h x, y), 

wherein at each coordinate (x, y), backscatter a 0 ((|) z , t z -, x, y) at 
a posting (x*, y^)=(0,0) is composed of the mean part a 0 (x, y), 
averaged over all (|) z and all t z , and a zero-mean fluctuation part 
of €(4> z -, t., x, y) variable in (|) z . at different time t z .. 


N 

G(x, y) = Yj 

i = 1 


G(4h,x, y) 
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